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Abstract 

Ion transport in biological and synthetic nanochannels is characterized by phenomena such as 
ion current fluctuations and rectification. Recently, it has been demonstrated that nanofabricated 
synthetic pores can mimic transport properties of biological ion channels [P. Yu. Apel, et al, 
Nucl. Instr. Meth. B 184, 337 (2001); Z. Siwy, et al, Europhys. Lett. 60, 349 (2002)]. Here, 
the ion current rectification is studied within a reduced ID Poisson-Nernst-Planck (PNP) model 
of synthetic nanopores. A conical channel of a few nm to a few hundred of nm in diameter, 
and of few ^m long is considered in the limit where the channel length considerably exceeds 
the Debye screening length. The rigid channel wall is assumed to be weakly charged. A one- 
dimensional reduction of the three-dimensional problem in terms of corresponding entropic effects 
is put forward. The ion transport is described by the non-equilibrium steady-state solution of the 
ID Poisson-Nernst-Planck system within a singular perturbation treatment. An analytic formula 
for the approximate rectification current in the lowest order perturbation theory is derived. A 
detailed comparison between numerical results and the singular perturbation theory is presented. 
The crucial importance of the asymmetry in the potential jumps at the pore ends on the rectification 
effect is demonstrated. This so constructed ID theory is shown to describe well the experimental 
data in the regime of small-to-moderate electric currents. 

PACS numbers: 05.60.Cd, 05.40.Jc, 81.07.De 
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I. INTRODUCTION 



The interaction of a cell with the extracellular environment through the membrane is of 
prominent importance for the cell behavior and cellular function. The lipid bilayer of cell 
membranes constitutes a barrier to the passage of charged and polar molecules [lj]. The 
channel proteins, which form narrow hydrophilic pores, primarily allow for the passage of 
small inorganic ions. The opening and closing of ion channels characteristically depend on 
the membrane voltage or the binding of selective ligand molecules j^j]. 

The production techniques and the application of single nanochannels in polymer films 
with lengths of micrometers and diameter sizes on the nanoscale is presently attracting 



widespread interest 



3j- Synthetic nanochannels have been fabricated as mimics of real 



biological nanopores and ion channels 4J. A recent intriguing finding was that conically 
shaped nanopores can be put to work as promising sensing elements for small molecules [5( , 
DNAs |6|, proteins [7| and yet other substances [8j. 

The sensing and transport properties of conical nanopores are strongly dependent on the 
shape of the pores, i.e. on parameters such as the cone angle and its length. An important 
feature of charged conical nanopores in comparison to the cylindrical ones is that the voltage 
drop caused by the ion current is centered around the narrow tip 

The flow of ions through a nanopore can either be driven by concentration gradients or 
an electric field, or both. It can approximately be modeled by means of an electro-diffusion 
equation. Within this description, the electric field inside the pore is governed in a self- 
consistent manner by the ion concentrations via the Poisson equation. At the boundaries, 
both the ion (bulk) concentrations and the electrical potentials are externally fixed. This 
may generate a steady current flow. With this work, our objective is the derivation of 
the general solution of such a non-equilibrium steady-state problem. A challenge is then 
the investigation of the influence of an externally applied voltage on different polarity for 
rectification, being induced by the asymmetry of the nanopore. 

A three-dimensional PNP modeling of the ion conductance in biological gramicidin A 
channels has been put forward in Ref. lo| for several membrane electrostatics. This 3D 
modeling was next investigated in Ref. [ll| in the context of the so called Potential-of- 
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Mean-Force-Poisson-Nernst-Planck theory which accounts for the dynamical relaxation of 
the channel forming protein and the surrounding medium by incorporating the free energy of 
inserting a single ion into the channel. Furthermore, OmpF and a-hemolysin porin channels 
were studied within similar approaches, complemented by molecular dynamics, and Brown- 
ian dynamics simulations in Ref. 



12, 



13] , and standard methods of continuum electrostatics 



recalled in Ref. Q. In part icn.hu:. for these porin channels, an asymmetric current-voltage 
characteristics, which implies a current rectification under alternating voltage conditions, 
has been observed. However, the structural complexity of biological ion channels and, in 
particular, an inhomogeneous charge distribution on the channel wall prevented any possibil- 
ity for a rigorous analytical treatment in these numerical studies. What is more, it remains 
unclear until now how a very popular one-dimensional, textbook 15j description of the ionic 



conductance of such and similar biological channels can be justified from 3D picture if the 
channel's length L and diameter D does not differ much, like L/D ~ 3 for porins. 

Synthetic nanochannels 3j, which have microns in length, are very much different in this 
respect providing a nice testbed for a one-dimensional reduced description. Such channels 
are wide enough (up to several nanometers in the narrowest part), so that the self-consistent 
mean field PNP theory, which neglects totally the correlation effects (existing e.g. due to a 
finite size of ions), seems well applicable. Moreover, the conductance of an ion channel in the 
imit, where its length exceeds largely the screening Debye length, should be strictly Ohmic 
1^ | . if the diameter does not vary and the channel wall is not charged. If the mobilities 
of free cations and anions are equal, the rectification effect can therefore only emerge due 
to a synergy of the entropic driving force caused by the channel asymmetry (conical pore) 
and the electrostatic effects caused by the fixed charges present on the channel wall. An 
analytical theory for such a system is currently absent. 

With this work we shall address this challenge within a singular perturbation theory 
approach for long conical pores, where the pore length substantially exceeds both its diameter 
and the Debye length. The obtained analytical results are validated by precise numerics. 
Moreover, numerical results are obtained for realistic surface charge densities beyond the 
validity range of perturbation theory. Then, the satisfying agreement with experimental 
data is demonstrated. 
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FIG. 1: (Color online) A sketch of the profile of the conical nanopore. The wall inside the pore 
exhibits a surface charge density a. The local radius of the pore, R(z), is given by Eq. ([5|). The 
drawing does not reflect the actual physical (z : r)-proportions. 

II. IONIC TRANSPORT 

Three dimensional Poisson-Nernst-Planck (PNP) equations 



We are interested in the ionic currents through relatively narrow pores which are created 
in a thin (of /im) dielectric film. The channel wall is assumed to be weakly charged by 
a prescribed surface charge density cr, (see Fig. [I]). The whole system is immersed in an 
electrolyte solution. We consider the limit of long channels with the length being much larger 
than both, the pore diameter and the Debye screening length: £d = \J e§e w k^T j (2e 2 A^4l0 3 / c ) 
(in meters). Therein, eo denotes the dielectric constant of vacuum, e w ~ 80 the relative 
dielectric constant of water, k& the Boltzmann constant, T the temperature, e the elementary 
charge, and Na is the Avogadro number. J c = J2i=± z^ 2 c;(oo)/2 defines the ionic strength 



of the electrolyte 



151 ] . It is given in terms of the bulk concentration of anions, c_(oo), and 



cations c + (oo) (in molars = moles/liter, M), respectively, with the valences being v\. For 
monovalent electrolytes i.e. V{ = ±1, with c + (oo) = c_(oo) = 1 M, the ionic strength is 
I c = 1 M and thus £d ~ 0.305 nm at room temperature T = 298 K. This presents a typical 
experimental situation for the synthetic pores studied in Ref. j], U, E], where £d ~ 1 nm at 
sub- molar concentrations of ca. 0.1 M. Our electrolyte system consists of two kinds of ions, 
potassium, K + , and chloride, Cl~. 
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Let us next start from the 3D electro-diffusion equation 

dci(r,t) - 

m =-V-Mr,t) (f) 

where 

ji(r, £) = -AVQ(f, £) + evifii£(r)Ci(r, t) (2) 

defines the mass flux density ji of the i-th species of ions, r is the position vector, Di the 
diffusion constant, \x\ the mobility of the ion particles, and 8 denotes the electric field. The 
consistency with the thermal equilibrium demands that the mobility of the ions /ij fulfills the 
Sutherland-Einstein relation, reading Di = fit/ (3, with the inverse temperature (3 = 1/ksT. 
Because £(r) = — V$(r), where $(r) is the electric potential, one can recast Eq. (CQ) into 
the more convenient form: 

dCl ^ l) = Di V exp[-ez/i (3 $(f)] V exp[ (3 $(r)] Cj(r, t). (3) 

The appropriate choice of coordinates for the problem at hand are the cylindrical coordinates 
(r, 0, z). The electric potential $(r, ^) is governed in a self-consistent manner by the Poisson 
equation 

e V • [e(r, z) V$(r, z)] = - Pi(r, z) - p fix (r, z). (4) 

i 

The ion concentration c, (in molars) is related to the density of electric charge pi (in units 
of Coulombs per cubic meter, C/m 3 ) through pi = 10 3 eui Na @ (—t + R(z)) Ci, where 6 is 
the Heaviside step function, and 

R(z) = R(0) + (R(L) - R(0))z/L (5) 

denotes a variable cone radius. The nanochannel radii at the left and right pore ends 
are R(0) and R(L), respectively, and L is the length of the pore. The charge density 
Pfix = b(j — R(z))a, with 8(x) being the Dirac delta-function, represents the fixed charges 
located on the inside of the channel wall. e(r, z) = e p 6(r — R(z)) + e w Q(— r+R(z)) describes 
dielectric properties of the open nanopore, where e p , e w are the relative dielectric constants of 
the polymer and water, respectively. The set of relations ([3]) and (jl]) for V-j'j = constitutes 
the system of coupled Poisson-Nernst-Planck (PNP) equations. 
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III. ONE-DIMENSIONAL MODEL REDUCTION FOR LONG PORES 



The primary purpose of our work is to gain analytical insight into the nature of rectifica- 
tion current. Since an analytical solution for the 3D PNP in general could not be found, a 
simplification is unavoidable. The ion diffusion and the ion flow are confined by the channel 
to a variable cross-section area nR 2 (z). Assuming local equilibrium in the transverse (r) 
direction, the full three dimensional system can be cast into an effective one dimensional 
problem. The price to be paid is the emergence of an additional, entropic potential contri- 



bution AS(z) 



17 
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20| | . Moreover, as pointed out by Zwanzig 17], a renormalization 



of the diffusion coefficient emerges as a dynamical effect due to a finite relaxation time in 
the transversal direction. In our situation (i.e. for a small opening angle of the cone) this 



latter dynamical correction to the diffusion coefficient is negligible 



2JJ . If the resulting cur- 



rents are too large these model assumptions can not be justified, as the particles must have 



20 1 



sufficient time to relax into the transverse direction 

In cylindrical coordinates (r, 0, z) the equation for mass current due to the coaxial sym- 
metry assumes the form 



d 

ji{r,z,t) = - Aexp[-ez/i/3$(r,2:)] — exp[ez/; /5$(r, z)\ c^r, z, t) r 

9 

- Di expf-ez/, (3 $(r, z)\ — exp[ez/i (3 $(r, z)\ Q(r, z, t) z, 

where f and z are the units vectors in cylindrical coordinates. 
The steady-state regime is described by the condition: 

V-jI = 0. 



(6) 



(7) 



Because the length of the channel by far exceeds its diameter, we can assume local 
equilibrium, i.e. uniform distributions for the ions and a constant potential across each 
channel cross-section, 



c(z, r, t) w c(z, t), $(r, z, t) w $(z, t). 



(8) 



The validity of an instant local equilibrium in transversal direction is fulfilled for weakly 
charged channel walls and not too large ion flow through the channel. Under this 



assumption, it is possible to reduce the complex original 3D problem to a simpler, effective 
one-dimensional description. 



Using the divergence of a given vector u 

V ■ u = lim - . . - — (j) u ■ h dS, (9) 

a 2 ^o A(z)Az Js 

where V(z) = A(z)Az denotes the volume surrounded by the closed surface S with the 
normal vector n, and the fact that no ion flux across the channel wall occurs, Eq. (TJ|) can 



be reduced 



221 ] under the assumptions in Eq. (JH]) to the form 

ddi(z,t) 8 f <9$(z) . , 8 Ci(z,t) 



9 [ n ^ , ,d<S>(z) A , , „ d Ci{z,t)} , . 

^PDrc^t^ + Mz)^-^, (io) 



dt 

wherein Ci(z,t) = 10 3 q(z, t)A(z) denote the one-dimensional concentrations (in moles per 
meter). 

Likewise, using (Q and the electrostatic boundary condition on the channel wall, reading 

e p£polymcr — = (H) 



the 3D Poisson equation is effectively reduced 



22 



23j into the form 



d 2 $(z) 2a , T cAzA) d$(z) d , , . , . N 

= -sw - ^ E ^ - ln (12) 

The quantity 2a/R(z) has the meaning of the volume density of the fixed surface charge, 
p = AQ/AV, where AQ = 2naR(z)Az and AV = nR 2 (z)Az. 
Likewise, Eq. (11 01) can be recast as 

8c (z t) 8 8 ( 8 1 

= -Q^Mz,t) = Tz I A exp[-0Gi(z)] ^SCM) eMPGM J, (13) 

where Jj denotes mass fluxes, and Gi(z) is the free energy (potential of mean force), given 
by 

G l (z) = -TS(z) + ev l <S>(z). (14) 
Due to the reduction from 3D into ID it inherits the entropic contribution 

S(z)=k B \n(A(z)/A ), (15) 



where Aq denotes an arbitrary, but irrelevant reference cross-section area. The coupled 
nonlinear equations ffT2l together with the set of Eqs. ffT3l - ffT5l) present a self-consistent, 
reduced ID electro-diffusion equation which is similar in spirit to the well-known reduction 



13, 



18 



19 



20|. The 



of a 3D-diffusion problem to an effective Fick- Jacobs diffusion equation 
dependence R(z) is still arbitrary within the discussed approximations; its variation must 
be smooth, however, such that the z-variations of the normal to the channel's surface from 
the cylindrical geometry become negligible up to the second order. This constructed ID 



Poisson-Nernst-Planck modeling generalizes the one given in Ref. 
pores to nanopores possessing a variable cross-section diameter. 



22j for strictly cylindrical 



A. Evaluation of the electric currents 

We consider the steady-state, non-equilibrium solutions characterized by constant mass 
fluxes Jj. The corresponding electric currents are ij = -FVjJj, where F denotes the Faraday 
constant. Explicitly, these electric currents for the i — th ion species read: 

h = -FDiUiexpl-pGiW—Ciiz) exp[(3G. t (z)] . (16) 

Thus, currents are fully determined by the bulk ion-concentrations and the electric potential 
difference U = $(0) — across the membrane, as determined by the effective lD-Poisson- 
Nernst-Planck equations. 



B. Dimensionless equations 

For the sake of convenience we perform all our numerical and analytical calculations in 
dimensionless variables. In doing so, we transform the Poisson-Nernst-Planck equations 
to their dimensionless form by use of following relations: z* = z/L, R*(z) = R(z)/R(0), 
= Pe$(z), I* = IiL/ (FDiCo), c*(z) = Cj(z)/c , where c is a reference ID-concen- 
tration (in 1/6.023 x 10~ 14 mole/m). Then, the equation for the constant electric currents 
assume the appealing form 

d . ,d<f>* . 2 dR* 
^ + ^ + ^—-,>-— = , (17) 



and (1T2]) transforms into: 

J_ (<P&_ 2_dR^d^\ _e_ 



with 



c F /?e L 2 

and the effective dimensionless surface charge density 

2aR(0) 



1/A 2 = €0€ ™T { rl (19) 



(20) 



c F 

Below, upon simplifying the notations, we shall suppress the notation with the super 
script (*). 



C. Boundary conditions 

We next solve the system of equations ( 1TTI) and ( fl8l) with the following boundary condi- 
tions: 

ck+(0) = VTC K +,L, Cd-(0) = 7TC C 1-,L, 

c K +(l) = vr(l + 7) 2 ck +)jR , cci-(l) =7r(l + 7 ) 2 c cl - ifl , (21) 

$(0) = 0, $(1) = $ R . 

wherein 7 = - R(z = 0))/R(z = 0) = i2(l) - 1, c h{L:R} = 10 3 c hulK{L , R} R 2 (z = 0)/c 

and Cbuik,{L,_R} denote the bulk concentrations of ions on the left and right sides, respectively. 
The electro-neutrality condition yields: ck+,l = cc\-,l — c l, and ck+,r = cci- t R = cr. 
Furthermore, the difference of dimensionless potentials across the nanopore is related to the 
applied voltage U (in units of Volts) by $(0) - $(1) = (3eU, yielding U = — <&(l)/(/3e). 

IV. SINGULAR PERTURBATION THEORY 

The corresponding ID system contains two small parameters, namely 

(i) 1/A, in Eq. (fl9l) . which is proportional to the ratio of the Debye length £d to the channel 
length L 

10 



(ii) the scaled surface charge density e [see Eq. (1201) and Appendix]. 

In our case i.e. a long nanopore with a typical length L = 12000 nm, 1/A is of the order 
of 10 -5 for sub-molar concentrations. The parameter e is of the order of 10 _1 , or larger. 
Therefore we will focus on the leading term in the limit 1/A — > 0, while e will be considered 
as a regular series expansion parameter. Thus, the problem requires to use a singular 
perturbation theory in 1/A and a regular expansion in e. We use the standard method of 



24| and seek for an approximate solution of Eq. f[T7|) and 



matched asymptotic expansions 
(JTSJ) in the form 

$ = $(o) + + . . . ; 

A 

c K+ = 4°1 + ±<£> + . . . , (22) 

(0) , 1 (1) 

cci- = c^L + jc" c >_ + .... 
Since the small parameter 1/A 2 appears in Eq. (|T8l) in front of the derivatives, we are dealing 
with a typical singularly perturbed boundary-value problem. 

Below we consider a symmetric bulk situation with equal ion concentrations on both mem- 
brane sides, cl = cr = c. After cumbersome calculations the uniformly valid approximation 
on [0, 1] is found in the first order of e to read: 



0>{z) « &°\z) = $0(2) -e^(Xz) + ^(z) + 0f (A(l -z)) + y c + ^iV^) } 



+ 0(e 2 ), (23) 



C{k+,ci-}(^) = ^c Ej0 (^) - j ± ir/R(z) + c SA (z) 

T 27rc0f (Az) =F 2ttc(1 + 7 ) 2 0f (A(l - z)) T n(2 + 1 )\+ 0(e 2 ), (24) 
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where the upper sign refers to K + ions, and the lower one to CI , respectively, and the 
expressions 

-Jo 



*o(z) 



jR(z) 
h 1 



- G, ) R 2 (z) = C R 2 (z) 



(25) 



7C0 R(z) 
irl 



+ E = <S> 



R- 



R(z) 



Ji 



2 7 C R*(z) jR(z) 



+ C 1 )R 2 (z) 



-7T$ R ^(Z - l)z, 



Tlln 



1 Ji 1 

+ 



+ 



C 2 7 \Q-fC R 3 (z) ' 2 7 R 2 (z) ' R{z) 
W7 +El , 



CojR(z) 

contain eight constants C , Ci, I , Ii, J , Ji, E , E x which are determined from the matching 
conditions (see in the Appendix). Furthermore, the left/right boundary layer potential 
variations in Eqs. fl23l) . ff24|) are given, respectively, by 

(26) 



<p[(\z) = — (exp[-V2c\z] - 1 



*f(A(l-*)) 



1 



exp[ 



^(1 -z)\-l 



2c(l + 7 ) 

where \/2c\/L = 1/£d, with £d being the Debye length in bulk. 

Assuming that the both sorts of ions have equal diffusion coefficients 
Dei- = D, the two constants J = —Ik+ + -^ci- an d I = Ik+ + ^ci- ( see Ref. |25j) have the 



26Li.e. D K , 



following physical meaning: J = J — eJ\ + 0(e 2 ) yields an approximation to the negative of 
the total mass current, whereas I = I — eli + 0(e 2 ) approximates the total electric current, 



/ % -2ttc(1 + 7 )$« ( 1 + t^T^—^R 



-24C1+7-; • (2?) 
Eq. (|27|) presents a main analytical result of this paper. In the original physical variables it 
reads 



x 



Ik+ + lex- « 27rFDc bulk (R(Q)R(L)/L) (eU/k B T) 
I a ( I I \ eU~ 



(28) 



where Cbuik is the concentration of ions in bulk. 
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It shows that the rectification effect appears already in the first order of e. The effect 
vanishes, when the pore is not charged, a = 0. This is because the rectification effects for 
the monovalent cations and anions for an uncharged pore are (trivially) "counter-directed" 
and compensate each other exactly if both sort of ions have equal diffusion coefficients. 
Obviously, for the cylindrical pore with 7 = 0, it goes away as well. Apparently, the current 
rectification is thus due to a synergy of the entropic effect caused by the pore asymmetry and 
the surface charge present. Moreover, one can attribute this result also to different volume 
charge densities 2o~/R(z) of the fixed charges at the both channels ends. Given this latter 
interpretation, one can expect also rectification effect for cylindrical pores inhomogeneously 
charged with two very different charge densities at the ends. 

Using the computer algebra system MAPLE, we obtained also analytical results for elec- 
tric currents and fluxes (total and individual) up to fourth order in e. The results are, 
however, rather cumbersome and of limited analytical insight, and thus are not explicitly 
displayed here. We evaluated them though numerically, see the corresponding figures given 
below. Furthermore, the PNP-system in Eqs. ffTTl) . (fl8l) with the boundary conditions in 
Ec, CD ™s integrated — >. y bv maklng use of a collocating .nethod with adaptive 
meshing |27|]. The results of the analytical perturbation theory and the numerical solutions 
compare very favorably where the perturbation theory is expected to work properly. 



V. RESULTS AND NUMERICAL COMPARISON 
A. Perturbation theory vs. numerics 

We consider channels with some fixed left opening radius R(z = 0) = 3 nm and two 
different right opening radii: 

- "long pore": L = 12000 nm , R(L) = 220 nm 

- "short pore": L = 200 nm , R(L) = 6.616 nm . 

In both cases the opening angle of the cone is identical, ip = 2 arctan[(i?(L) — R(z = 0))/L] w 
2°. Furthermore, two different surface charge densities are used in our studies, namely: 
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a = —0.02, —0.1 e/nm 2 . The first one corresponds to a parameter value e = 0.12 which is 
well within the regime of validity of the perturbation scheme. For the second value we have 
e = 0.6 for which we already expect the perturbation method to fail, but it still might work 
occasionally. Furthermore, the diffusion coefficients were taken as = Dq\- = 2 x 10 9 
nm 2 /s, see in 26j, temperature T = 298 K and e w = 80. The concentration profiles for 

0.22 

OIK , , , , A I A 

a) 




0.22 




100 
z [nm] 



0.1008 




0.1004 - 



11990 11995 , 12000 



4000 8000 
z [nm] 



12000 



FIG. 2: (Color online) Concentration profile c K +(z) at equilibrium (c^ = cr = 0.1 M, $(z = 0) = 
$(z = L) = 0). Calculations are done for the "short pore" (a) and the "long pore" (b). Solid line 
and squares: a = —0.02 e/nm 2 , dashed line and circles: a = —0.1 e/nm 2 . Symbols in all cases 
denote the results of the numerical solution, lines represent the results of the perturbation theory. 
The insets depict a closer look into the left and the right boundary layers, respectively. 



both, the singular perturbation solution and the numerical solution are depicted in Fig. [3 



14 



The first order approximation agrees very well with numerical solution for the small surface 
charge density of o = —0.02 e/nm 2 . However, this agreement is worsening upon increasing 
the values of a. For the moderately large charge density o = —0.1 e/nm 2 the discrepancy 
between the exact solution and analytical approximation becomes already significant. In 
the inset in Fig. [2] we show the left and the right boundary layers, according to Eq. (1241) . 
Both layers possess a width of several Debye lengths In the case of the "long pore", 
the increase of the concentration ck+{z) within the boundary layer is more distinct near 
the narrow opening (the boundary z = 0), as compared to the wide opening at the right 
boundary at z = L, see Fig. [2b. This is so because for this channel the right opening radius 
R(L) = 220 nm is very wide. The shorter channel with left and right opening radii of similar 
size does not display such a striking difference. 

B. Current- Voltage characteristics and rectification 

In Figs. [3] and H] we depict the current- voltage (I-U)-characteristics of the ion transport. 
The comparison of analytical and numerical results show that the currents for both the 
"short pore" (Fig. [3h) and the "long pore" (Fig. [3b) are well predicted by the first order 
perturbation theory, cf. Eq. (I2"gj) . in the range of U from —0.1 V to 0.1 V, for both surface 
charge densities a = —0.02 e/nm 2 and o = —0.1 e/nm 2 . However, for the larger density, 
the agreement deteriorates for the "short pore" . As expected, the discrepancy grows further 
with increasing the absolute value of the applied voltage U across the nanopore, cf. Fig. [H 

The above examples evidence that the first order perturbation expansion works well only 
for a relatively small charge density o and not too large applied voltages. However, in a wider 
range of voltages and for larger charge densities one has to use a higher order perturbation 
theory in e. 

In Fig. O we depict the numerical I — U current-voltage characteristics for both the 
potassium Jk+, and the chloride Iqi- currents versus the fourth order perturbation theory 
results. The found agreement is rather good. Undoubtedly, the characteristics are strongly 
non-linear and asymmetric. Consequently the nanopore exhibits rectification properties. We 
can analytically describe this effect in the first order of e. Then, the individual ionic currents 
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FIG. 3: (Color online) Current voltage (I — U) characteristics of the "short pore" (a) and the "long 
pore" (b) for U from —0.1 V to 0.1 V. Solid line and squares: a = —0.02 e/nm 2 , dashed line and 
circles: a = —0.1 e/nm 2 . Symbols denote the numerical solution, the lines represent the results of 
first order perturbation theory (in e). 

are approximated by the following expansions: I K + = (I — Jo)/ 2 — e(Ii — J\)/2 + 0(e 2 ) 
and J C1 - = (J Q + Jo)/ 2 — e(ii + J\)/2 + 0(e 2 ). Interestingly, for both the potassium and the 
chloride currents the difference between the absolute values of positive and negative current 
branches is the same and equals to: 

|/ {K+) ci- } (C/)| - |/ {K+ ,ci-}(-C/)| - -e/i = - l -D.a R{L) ~ m . (29) 

In a negatively charged nanopore, cations experience a potential well, whereas anions sense 
a potential barrier. Thus, the potassium concentration and also the potassium current 
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FIG. 4: (Color online) Current voltage (I — U) characteristics of the "short pore" (a) and the "long 
pore" (b) for U from —0.5 V to 1 V. Solid line and squares: a = —0.02 e/nm 2 , dashed line and 
circles: a = —0.1 e/nm 2 . Symbols in all cases stand for numerical solution, lines represent first 
order perturbation theory (in e). 

are always larger than those of chloride, see Fig. This effect becomes enhanced upon 
increasing the absolute value of the (negative) charge density a. However, the equivalent 
asymmetry of the I — U dependence for individual currents indicates that the rectification 
quality of the pore is independent of the ion sign within this order of perturbation theory. 

Furthermore, in Fig. E^i we compare the electric potential profiles $(-z) of the "short pore" 
(squares) (L = 200 nm) and the "long pore" (L = 12000 nm) (circles) for the equilibrium 
situation. They coincide up to z = 195 nm, because the angle of the cones of both nanopores 
are the same. This implies that the electric field profiles in the narrow parts of both channels 
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FIG. 5: (Color online) The potassium and chloride currents in the "short pore" for a = —0.02 
e/nm 2 . The solid line represents the fourth order (in e) perturbation calculation for potassium 
and the dashed line for the chloride current. Symbols in all cases present the numerical solution 
(circles: potassium current; squares: chloride) 



are identical. Since the rectification efficiency seemingly is caused mostly by the sharp 
potential variations in the narrow regions, one might naively expect that the rectification 
quality will be similar for both cases, even if the absolute I — U characteristics for these two 
nanopores are rather different, cf. Fig. (6h>. However, this would be an incorrect conclusion. 

Indeed, the absolute values of the total electric current I = I K + + I C \- is lower for the 
"long pore" as compared to the "short pore". This result corroborates with the higher 
resistance of the longer pore. To quantify the rectification effect, we define the quantity 

° = |/(£0l + |/(-£0l <30) 



as a rectification-measure. For it, the lowest order perturbation theory result in Eq. (l28j) 
yields 

_ I I ( 2a 2_^_\e\Ul 

24c bulk F \R(0) R(L)J k B T' ° ' 



As seen from this expression, the quality of rectification a, does not depend on the chan- 
nel length L and is determined rather by the difference of the inverse radii and by the 
surface charge. Note that, in the above-described model, the volume density is given by 
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a(z) = 2a/R(z) (Eq. [12]), where R(z) denotes a variable cone radius. The geometry of the 
charged conical pore results not only in the presence of entropic potential in the ID current 
equation (Eq. [TBI) but also in the asymmetrical volume charge distribution a(z) in the re- 
duced ID Poisson equation (Eq. [T2l . This fact gives strong indication that in conical pores 
the rectification is mainly due to different potential jumps in the boundary layers at both 
ends of the channel. This is in line with the previous statement that the I — U characteris- 
tics of the nanopore depends crucially on the total asymmetry of the potential profile [28]. 



Accordingly, a reveals, cf. Fig. [6b, that the "long pore" acts as a better rectifier, despite its 
larger resistance. However, this is not due to a larger length, but due to a larger R(L). 



Theory vs. experiment 

r~i 

In Fig. [3, a comparison between the experimental data [29J] and theoretical results is 
presented for two different nanopores of length L = 12000 nm with the radii -R(O) = 3.5 
nm, R(L) = 340 nm and -R(O) = 3 nm, R(L) = 265 nm, respectively. The experimental 
charge density is not known precisely. The assumed value is around a = —1.0 e/nm 2 , 
which corresponds to e ~ 6. This parameter value is clearly beyond the validity range of 
perturbation theory. Nevertheless, one can use numerics. Given the approximate, reduced 
character of the studied lD-Poisson-Nernst-Planck model, the semi-quantitative agreement 
between the theory and the experiment is quite satisfactory. 

Indeed, the agreement is rather satisfactory for not too large currents, / < 0.5 nA. 
The simplest explanation for the discrepancy occurring for strong currents could be the 
violation of the condition of a local equilibrium in the transverse direction (see Sec. One- 
dimensional model reduction). Then, this breakdown of the ID description would be similar 
in nature to the one observed for diffusive transport of biased, non-interacting particles 



m 
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211 ] . To fully justify this statement, one can follow the analysis in Ref. [20] to find 
the different time scales characterizing the problem, i.e. the time scale for diffusion in the 
transverse direction and the time scale associated with the drift of the i-th ion along the axial 
direction. By comparison of this two time scales one could find out whether transported 
ions have enough time to relax in the transversal direction. However, in the present, more 
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FIG. 6: (Color online) In the panel (a) we depict the potential profile $(z) at equilibrium (cl = 
cr = 0.1 M, $(0) = $(L) = 0) for a = -0.02 e/nm 2 . In the panel (b) we depict the I - U 
dependence and in panel (c) we show the rectification measure a. The calculations are done for 
the "long pore": dashed line and circles; and the "short pore": solid line and squares. The symbols 
in all cases stand for numerical solution; the lines represent the perturbation theory (in e): up to 
second order (a), and in (b) and (c) up to fourth order. 

complicated situation the proposed consideration results in a very rough and somewhat 
vague estimation (because of a heterogeneous distribution of the ions and the corresponding 
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FIG. 7: (Color online) The current voltage (I—U) dependence in the experiment (symbols) and the 
1D-PNP theory (lines). Only the numerical results are given because the perturbation approach 
fails for the large surface charges used in the experiment. We depict results for one value of the 
surface charge density a = —1.0 e/nm 2 and for cl = cr = 0.1 M which lies within the range of 
experimentally measured values. We use two different pore geometries: panel (a) R{z = 0) = 3 
nm, R(L) = 265 nm, L = 12000 nm; panel (b) R(z = 0) = 3.5 nm, R(L) = 340 nm, L = 12000 
nm. 

electrical field along the pore). For better agreement a 3-D PNP modeling would be called 
for. 

Such a 3D PNP study, would in natural way, reveal the dependence of the electric po- 
tential on the radial coordinate r, and, in turn, a heterogeneous distribution of ions in the 
transverse direction, contrary to the model assumptions in Eq. (JHJ). Clearly, the larger is 
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the surface charge density, the stronger will be the violation of these model assumptions. 
Because of a negatively charged channel wall, the positive ions should adhere to it forming 
an electrical double-layer of the width £d with an exponentially enhanced concentration of 
cations. A physical criterion to disregard this complexity can be obtained by demanding that 
the electrostatic energy per ion near the channel wall (in the absence of external voltage), 
eV ~ ecx^D/(eoe«)), is smaller than k^T. For e w = 80 and £d = 1 nm this yields \a\ < 0.1 
e/nm 2 . For higher surface charge densities (like in the experiment) not only the homogeneity 
assumptions become increasingly violated, but also the surface current contributions 30l. |31| 



cannot be neglected. These effects are beyond our description within a 1D-PNP reduction. 



15, 
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Moreover, for higher concentrations c ~ 1 M the ion-size effects become important 
implying the presence of ion correlations. This goes beyond the mean field description. 

To summarize, a numerical study of the full three-dimensional set-up would be helpful in 
order to distinguish between the principal shortcomings of the PNP approach and the failure 
of our ID reduction for large currents and/or high surface charge densities. In particular, the 
overall failure of the PNP approach can be due to formation of the Stern layer of condensed 

n 

counter-ions on the pore wall [15]. Such a layer will screen and renormalize (i.e. decrease) the 
surface charge density "seen" by other mobile ions, thus reducing the "volume" rectification 
effect. A proper treatment of the surface current rectification effects in such a layer would 
bring about the theme of strong correlations in the ionic transport which is beyond the scope 
of this work. 



VI. SUMMARY AND CONCLUSION 

In this paper we provided an analytical treatment of the problem of current rectification 
by artificial conical nanopores made in synthetic membranes of about submicron-to-micron 
width. Within a one-dimensional reduction of the original three-dimensional problem, we 
have derived 1D-PNP equations, cf. Eqs. (12)-(14), which incorporate both the entropic 
effects and the charge density renormalization due to the variable diameter of the tube. 



These equations generalize the 1D-PNP equations in Ref. 
variable diameter. 



221 ] to the case of pores of a 
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Furthermore, we provided a singular perturbation treatment of the problem of ion con- 
ductance and rectification within the reduced 1D-PNP description implementing rigorous 
boundary conditions (12T|) . Our theory applies to nanopores with lengths exceeding largely 
the Debye length and the pore diameter. The developed theory corresponds precisely to the 
experiments done by Siwy et al. 29]. The validity range of perturbation theory requires, 



however, that the channel wall is charged weakly. Then, it agrees well with the numerical 
solution of the 1D-PNP problem for charge densities up to |cr| = 0.02 e/nm 2 . Unfortunately, 
the experimental charge densities are much larger, about \a\ ~ 1 e/nm 2 . Here the pertur- 
bation theory fails. However, the 1D-PNP equations can be integrated numerically and the 
obtained numerical solutions provide a good agreement with the experimental data by Siwy 
et al. 



29j for sufficiently small currents. 
Moreover, we quantify the rectification properties of conical nanopores. In the lowest 
order perturbation theory we obtained analytic formulas for the rectification-current / and 
the rectification measure a, see Eq. (130]) . The latter measure clearly indicates that the 
rectification property is caused by the difference in the volume charge density of the fixed 
charges. In other words, the rectification is due to an asymmetry in the potential jumps at 
the channel ends. 

A discrepancy between the experiment and the 1D-PNP modeling for large positive volt- 
ages most likely indicates the violation of the condition of local equilibrium in the transverse 
direction at strong currents, cf. Fig. [7J For this reason, the description within the 1D-PNP 
modeling breaks down and a full 3D-PNP treatment becomes necessary. We conjecture that 
the numerical solution of the 3D-PNP description will provide satisfactory agreement with 
the experimental data for the electrolytic solutions of mono- valence ions. 
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APPENDIX: STUDY CASE OF SINGULAR PERTURBATION THEORY 

We consider a set of two coupled differential equations (|17[) . ffTSj) with the boundary 
conditions given in Eq. (12 f I) . Each of these equations one can be presented in the form 

Px/xAv) = (A.l) 

where y = {ck+, cq-, and where two small parameters 1/A and e determine the behavior 
of the solution. The solution shall be denoted by y^ 1 ^'- The limiting problem (1/A = 0) 

PoM = (A.2) 

possesses the solution yi?\ The character of the problem changes discontinuously at 1/A = 0. 
This implies that we have a singular perturbation in this very parameter. 



The Outer Approximation 

In the outer region (corresponding to the interior of the 1-D channel) we approximate 
■ " ' ' ' — ^ c^ + (z), 5q_ (z), $(°)(z) j by use of the regular perturbation expansion: 

■ (0) ^ " Vo(z) + ey x {z) + e 2 y 2 {z) + . . . . (A.3) 



The Boundary Layers 

The region near the boundaries at z = and z = 1 wherein y is changing rapidly presents 
the boundary layer (left and right, respectively). The outer expansion looses validity there. 
We re-scale the problem near z = (z = 1) by setting ( = \z (x = X(z — 1)), and 
express the various functions |c^|(2;), (z), $^(z)| in terms of the new coordinates 
as {^(C; l/X),n L (C 1/A),0 L (C; 1/A)}, and {p R ( X ; l/X),n R ( X ; l/X)A R (x; 1/A)}, in the left 
and right boundary layers, respectively. Next, considering the limit expansions obtained by 
holding ( and x fixed and letting 1/A — > 0, the problem can be solved in terms of the regular 
expansion yf\e, z) and y^\e, z). 
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Matching procedure 

We choose four constants which allow that yz,(e, z) and y ou tcr(e, z), yn{e, z) and y ou tcr(e, 2), 
coincide for each order of the expansion in powers of e (as 1/A — > 0) in some intermediate 
zone between the left boundary layer and the outer region, and the right boundary layer 
and the outer region, respectively. This yields: 



Jo = 0, Co = 2nc, 

Jo = -C (l + 7 )$fl, E = Q-^Qr. 

7 

Ji= ^(2 + 7)^, C 1 = ^ R , 

To obtain an approximate solution y u that is valid uniformly on [0, 1] we finally add the 
boundary and outer approximations and subtract their common limit in the intermediate 
zone. 
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